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ABSTRACT 

Numerical simulations based on the ACDM model of cosmology predict a large 
number of as yet unobserved Galactic dark matter satellites. We report the results 
of a Large Area Telescope (LAT) search for these satellites via the 7 -ray emission 
expected from the annihilation of weakly interacting massive particle (WIMP) dark 
matter. Some dark matter satellites are expected to have hard 7 -ray spectra, finite 
angular extents, and a lack of counterparts at other wavelengths. We sought to 
identify LAT sources with these characteristics, focusing on 7 -ray spectra consistent 
with WIMP annihilation through the bb channel. We found no viable dark matter 
satellite candidates using one year of data, and we present a framework for interpreting 
this result in the context of numerical simulations to constrain the velocity-averaged 
annihilation cross section for a conventional 100 GeV WIMP annihilating through the 
bb channel. 

Keywords: dark matter - galaxies: dwarf - gamma rays: galaxies 


1. Introduction 

Astronomical evidence and theoretical arguments suggest the existence of non-baryonic cold 
dark matter (CDM) . In standard model cosmology, dark matter constitutes approximately 85% of 
the matter density and nearly one-fourth of the total energy density of the universe (Komatsu et al. 
2011). While very little is known about dark matter beyond its gravitational interactions, a popular- 
candidate is a weakly interacting massive particle (WIMP) (Jungman et al. 1996; Bergstrom 2000; 
Bertone et al. 2005). From an initial equilibrium state in the hot, dense phase of the early universe, 
WIMPs freeze-out with a significant relic abundance to constitute much, if not all, of the dark 
matter in the universe. In regions of high dark matter density, WIMPs may continue to annihilate 
into Standard Model particles through the same processes that originally set their relic abundance. 

Gamma rays produced in the final state of WIMP annihilation, either mono-energetically from 
direct annihilation or as a continuum of energies through annihilation into intermediate states, may 
be detectable by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope 
(Fermi). These 7 rays would be produced preferentially in regions of high dark matter density 
within the Milky Way or extragalactie sources. The LAT has reported upper limits on mono- 
energetic 7 -ray line from annihilations in the smoothly distributed dark matter halo surrounding 
the Galactic plane (Abdo et al. 2010c). Additionally, upper limits on the continuum 7 -ray emission 
from WIMP annihilation have been placed using dwarf spheroidal galaxies (Abdo et al. 2010e; 
Ackermann et al. 2011a), the diffuse isotropic background (Abdo et al. 2010a), and clusters of 
galaxies (Ackermann et al. 2010). 

Cosmological A r -body simulations predict that the Milky Way should have many more satellites 
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than are currently observed at optical wavelengths (Diemand et al. 2007; Springel et al. 2008). This 
prediction allows for the possibility that the majority of these satellites are composed solely of dark 
matter. While these simulations have the resolution to recover all satellites down to a mass of ~ 10 6 
Mq, the minimum mass of a bound satellite orbiting the Milky Way may be as low as the Earth 
mass (~ 10~ b Mq), a scale that is roughly set by the WIMP velocity at freeze-out (Hofmann et al. 
2001; Loeb & Zaldarriaga 2005). Massive dark matter satellites located in the outer regions of the 
Galactic halo or lower mass satellites located near to Earth may constitute significantly extended 7 - 
ray sources that could be detectable with the LAT. It has been suggested (Tasitsiomi & Olinto 2002; 
Koushiappas et al. 2004; Pieri et al. 2008; Baltz et al. 2008; Springel et al. 2008; Anderson et al. 
2010) that a significant 7 -ray flux from WIMP annihilation could arise from these dark matter 
satellites within the Milky Way halo. 

Here, we report on a search for dark matter satellites via 7 -ray emission from WIMP anni- 
hilation. We begin by examining the theoretical motivation for such a search. Then, we describe 
the selection of unassociated, high Galactic latitude 7 -ray sources from both the First LAT Source 
Catalog ( 1 FGL) (Abdo et al. 2010b) and an independent list of source candidates created with 
looser assumptions on the source spectrum. The likelihood ratio test was used to distinguish ex- 
tended sources from point sources and WIMP annihilation spectra from conventional power-law 
spectra. No dark matter satellite candidates were found in either the unassociated 1 FGL sources 
or the additional list of candidate sources. Finally, Via Lactea II (Diemand et al. 2007) and Aquar- 
ius (Springel et al. 2008) simulations were used to derive upper limits on the annihilation cross 
section for a 100 GeV WIMP annihilating through the bb channel. 


2. Dark Matter Satellites 

In this section, we describe models for the spatial and mass distributions of dark matter 
satellites in the Galactic dark matter halo, as well as models for the internal density profiles of 
these satellites. We motivate the search for dark matter satellites lacking optical counterparts 
by predicting the number of satellites that are detectable, given the constraints on dark matter 
annihilation from dwarf spheroidal galaxies. We discuss how extrapolating the mass function of 
dark matter satellites below the mass resolution of current numerical simulations impacts LAT 
detection potential. 


2.1. Numerical Simulations 

The Aquarius (Springel et al. 2008) and Via Lactea II (VL-II) (Diemand et al. 2007) projects 
are currently the highest resolution numerical simulations of dark matter substructure at the scale 
of Galactic halos. As part of the Aquarius project, six independent simulations of a Milky Way- 
mass dark matter halo were generated. The VL-II project provides an additional, independent 
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simulation. These simulations model the formation and evolution of a Milky Way-size dark matter 
halo and its satellites having over one billion ~ 10 3 M® particles. Each simulation resolves over 
50,000 satellites within its respective virial radius, which is defined as the radius enclosing an average 
density 200 times the cosmological mean matter density. Each bound satellite has associated with 
it a position with respect to the main halo, a velocity, a tidal mass, M t jdai , a maximum circular 
velocity, V^n ax , and a radius of maximum circular velocity, Rv mgx - 


Generally, the 7 -ray flux from annihilating dark matter in a satellite (0wimp) can be expressed 
as a product of the line-of-sight integral of the dark matter distribution in a satellite ( J-factor) and 
a component depending on the particle physics model (<F PP ) for WIMP annihilation (Strigari et al. 
2007) 

0wimp (E) = Jx $ pp (f?) ( 1 ) 

where 
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Equation (2) represents the line-of-sight integral along l through the satellite dark matter density 
profile p, where D is the distance to the satellite center, 6 is the offset angle relative to the center, 
and r{9, D, l ) = VD 2 + l 2 — 2 Dl cos 9 is the distance from the satellite center. The solid angle 
integral is performed over A ft = 2^(1 — cos0). For satellites at large distances from the Earth, the 
J-factor can be approximated by Equation (3) where the volume integration is performed out to 
the satellite tidal radius (Tyler 2002). In Equation (4), (av) is the velocity-averaged annihilation 
cross section, Mwimp is the mass of the dark matter particle, and the sum runs over all possible 
pair annihilation final states with dNf/d,E and Bj representing the photon spectrum and branching 
ratio, respectively. 


For the present analysis, it was important to model the dark matter distribution within the 
satellites themselves. A Navarro— Frenk— White (NFW) profile (Navarro et al. 1997) with scale 
radius r s and scale density p s was used to approximate the dark matter distribution in the satellites: 


p{r) = 


Ps 

(r/r s )( 1 + r/r s ) 2 


An NFW profile was uniquely defined for each satellite from the values of F m » v and Ry mix using 



the relations (Kuhlen et al. 2008) 
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where r s is in kpc and p s is in Mq kpc“ 3 . The existence of dark matter substructure within the 
satellites themselves would increase the J-factor, but this contribution is expected to be no greater 
than a factor of a few (Martinez et al. 2009; Anderson et al. 2010). Thus, we took a conservative 
approach and did not include this enhancement when calculating J-factors. 

Although the NFW model is widely used for its simplicity and is broadly consistent with the 
dark matter distribution in satellites (Springel et al. 2008), a few caveats are important. First, 
simulations with increased resolution have revealed that more scatter exists in the central dark 
matter density profile than is implied by the NFW profile (Navarro et al. 2010). In fact, these 
simulations indicate that the central dark matter density in satellites is systematically shallower 
than the r _1 central density implied by the NFW model. However, these deviations occur at a 
scale of < 10~ 3 of the halo virial radius (Navarro et al. 2010), and thus do not strongly affect the 
predicted flux given the LAT’s angular resolution. Second, satellites are, in nearly all cases, more 
severely tidally truncated than the r -3 outer density scaling of the NFW profile. For satellites with 
large r s , using the NFW profile will result in a slight (< 10%) overestimation of the predicted flux. 


2.2. Extrapolation to Low-mass Satellites 

The satellite mass functions for the Aquarius and VL-II simulations are complete down to 
~ 10 6 Mq. However, theoretical arguments suggest that the mass function of satellites may extend 
well beyond this resolution limit, perhaps down to Earth-mass scale ( 10 6 Mq). Therefore, it is 

likely that the satellites resolved in the Aquarius and VL-II simulations are only a small fraction 
of the total number of bound dark matter satellites present in our Galaxy. We refer to these 
unresolved dark matter satellites with M t jdal < 10 6 M© as “low-mass” satellites. 

To estimate the contribution of low-mass satellites to the LAT search, we extrapolated the 
distribution of satellites in VL-II down to 1 Mq. Assuming a power-law mass function for satel- 
lites, dN/dMtidal MjMj 50 (Madau et al. 2008; Springel et al. 2008), we calculated the number of 
satellites at a given Mtida! within 50 kpc of the Galactic center. These low-mass satellites were 
distributed within this 50 kpc radius in accordance to the radial distribution described in Madau 
et al. (2008). The cut at 50 kpc is conservative based on the null-detection of the Segue 1 dw r arf 
spheroidal with a tidal mass of ~ 10' M© at ~ 28 kpc from the Galactic center (Geha et al. 2009). 

In order to model the internal dark matter distribution of low-mass satellites with an NFW 
profile (using Equation (6) and (7)), we fit relationships between Af tidal, L,y . and Ry m . v . For 
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VL-II satellites within 50 kpe of the Galactic center and with Mtidal > 10 6 M©, we found that 

Y T r /-^tidalV 

with Vo = io _1 20±0 - 05 kms -1 , /3 = 0.30 ±0.01, and a log-Gaussian scatter of oy max = 0.063 kms -1 . 
Additionally, we found that 

= Sc (9) 

with Ro = 10 -3 ' 1±0 ' 4 kpc, S = 0.39 ±0.02, and a log-Gaussian scatter of VR Via!ix — 0.136 kpc. Using 
these relationships, we randomly generated low-mass satellites consistent with the VL-II simulation 
down to a tidal mass of 1 M©. 


2.3. Comparison to Optical Satellites 


The VL-II simulation, extrapolated as described above, provides a theoretical model for the 
population of Milky Way dark matter satellites from 1O 1O M0 to 1 M©. A simple estimate of the 
detectable fraction of these dark matter satellites can be obtained from the 11-month limits on the 
WIMP annihilation flux from dwarf spheroidals (Abdo et al. 2010e). No 7-ray signal was detected 
and the strongest limits on the annihilation cross section result from the analysis of the Draco 
dwarf spheroidal, which has a J-factor (integrated over the solid angle of a cone with radius 0.5 °) of 
~ 10 19 GeV 2 cm“ 5 . The central mass of Draco is well known from stellar kinematics (Strigari et al. 
2008; Walker et al. 2009; Wolf et al. 2010); however, the total dark matter mass for Draco is less 
certain due to the lack of kinematic measurements in the outer regions of the halo. A conservative 
lower bound on the total dark matter halo mass of Draco was taken to be 10 s M©. Based on this 
lower bound, we determined what fraction of the satellites have a larger detection potential than 
Draco, i.e., what fraction of satellites have a J-factor greater than that of Draco. 


For an NFW profile, the J-factor is roughly described by 


J oc 


9 

r sP% 


D 2 


oc 


M °- 81 
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( 10 ) 


using the relationships that r s oc M 0,39 (from Equations (6) and (9)) and p s oc M -0 ' 18 (from Equa- 
tions (6) — (8)). Equation (10) makes it possible to compare the relative astrophysical contribution 
to the 7- ray flux for different halos based on their tidal mass and distance from Earth (Figure 1). 
The choice of particle physics annihilation model merely scales all satellites by the same constant 
factor. 


Figure 1 serves as a guide for evaluating the detectability of low-mass satellites. While the total 
number of satellites increases with decreasing mass, the J-factors of these low-mass satellites tend 
to decrease. This means that low-mass satellites, while dominating the local volume in number, 
are a subdominant contributor to the 7-rav flux at the Earth. Using the procedure discussed in 
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Section 5, we verify that extending the VL-II mass function to low mass has a minimal effect (< 5%) 
when setting upper limits on (av ) ; consequently, we do not consider these satellites in our primary 
analysis. For low-mass satellites to dominate the y-ray signal, a mechanism must be invoked to 
either increase the concentration for low-mass satellites, or decrease the slope of the mass function. 
Of course, the above statements do not preclude the possibility that there could be a low-mass 
satellite with a high J-factor very near to the Sun. 

In the context of the CDM theory, several dark satellite galaxies with no associated optical 
emission could be detectable by the LAT. In addition to motivating our satellite search, Figure 1 
allowed us to narrow our focus to those satellites with the best prospects for detection. Using 
Equation (10). we omitted satellites with J-factors more than an order of magnitude less than the 
lower bound on the J-factor for Draco. This greatly reduced the number of satellites for which the 
full line-of-sight integral in Equation (2) was calculated. 

3. Methods 

In this section, we review the tests applied to LAT sources lacking associations in other wave- 
lengths to determine if any are consistent with dark matter satellites. First, we summarize our data 
set and give an overview of an independent search for LAT sources without spectral assumptions. 
Then, we define a procedure for selecting candidate dark matter satellites using the likelihood ratio 
test to evaluate the spatial extension and spectral shape of each source. The ability of the LAT to 
detect spatial extension and spectral shape depends on source flux and spectral hardness. Exten- 
sive Monte Carlo simulations were required to determine cuts for rejecting both point-like sources 
and sources with power-law spectra at 99% confidence. When combined, these cuts allowed us to 
select for extended, non-power-law sources in our sample of high-latitude unassociated LAT source 
candidates with a contamination level of 1 in 10 4 . The work presented in this section is largely 
derived from Wang (2011). 


3.1. Data Selection 

Fermi has been operating in sky-scanning survey mode since early 2008 August. The primary 
instrument on board Fermi is the LAT, designed to be sensitive to y-rays in the range from 20 MeV 
to > 300 GeV. The LAT has unprecedented angular resolution and sensitivity in this energy range, 
making it an excellent instrument for detecting new y-ray sources (Atwood et al. 2009). 

Our data sample consisted of “Diffuse” class events from the first year of LAT data collection 
(2008 August 8 to 2009 August 7), which overlapped substantially with the data used for the 
1FGL (Abdo et al. 2010b). To reduce y-ray contamination from the bright limb of the Earth, we 
rejected events with zenith angles larger than 105 ° and events taken during time periods when the 
rocking angle of the LAT was greater than 47 ° (the nominal LAT rocking angle was 35 ° during this 
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time period). Due to calibration uncertainties at low energy and the current statistical limitations 
in the study of the instrument response functions (IRFs) above 300 GeV, we accepted only photons 
with energies between 200 MeV and 300 GeV. This analysis was limited to sources with Galactic 
latitudes greater than 20°, since the Galactic diffuse emission complicates source detection and the 
analysis of spatial extension at lower Galactic latitudes. We modeled the diffuse 7 -ray emission 
with standard Galactic (glLiem-v02.fit) and isotropic ( isotropic-.iem~v02.txt ) background models 1 . 
Throughout this analysis, we used the LAT ScienceTools 2 version v9rl8pl and the P6.V3JDIFFUSE 
IRFs 3 . 


3.2. Source Selection 

The 1FGL is a collection of high-energy 7 -ray sources detected by the LAT during the first 
11 months of data taking (Abdo et al. 2010b). It contains 1451 sources, of which 806 are at high 
Galactic latitude (| 6 | > 20°). Of these high-latitude 1FGL sources, 231 are unassociated with 
sources at other wavelengths and constitute the majority of the sources tested for consistency 
with the dark matter satellite hypothesis. However, the 1FGL spectral analysis, including the 
threshold for source acceptance, assumed that sources were point-like with power-law spectra. This 
decreased the sensitivity of the 1FGL to both spatially extended and non-power-law sources, which 
are characteristics expected for dark matter satellites. In an attempt to mitigate these biases, we 
augmented the unassociated sources in the 1FGL with an independent search of the high-latitude 
sky. 

We performed a search for 7 -ray sources using the internal LAT Collaboration software pack- 
age, Sourcelike (Abdo et al. 2010d; Grondin et al. 2011). Sourcelike performs a fully binned 
likelihood fit in two dimensions of space and one dimension of energy. When fitting the spectrum 
of a source, Sourcelike fits the fraction of counts associated to the source in each energy bin 
independently. The overall likelihood is the product of the likelihoods in each bin. This likelihood 
calculation has more degrees of freedom than that performed by the LAT ScienceTool, gtlike, 
which calculates the likelihood from all energy bins simultaneously according to a user-supplied 
spectral model. In this analysis, we used 11 energy bins logarithmically spaced from 200 MeV to 
300 GeV. 

Using Sourcelike, we searched for sources in 2496 regions of interest (ROIs) of dimension 
10° x 10° centered on HEALPix (Gorski et al. 2005) pixels obtained from an order four tessellation 
of the high-latitude sky (; 6 | > 20°). Each ROI was sub-divided into 0.1° x 0.1 0 pixels, and for each 
pixel the likelihood of a point source at that location was evaluated by comparing the maximum 


1 http:/ /fermi.gsfc. nasa.gov/ssc/data/access/lat/BackgroundModels.html 

2 http://fermi.gsfe.nasa.gov/ssc/data/analysis/software 

3 http:// fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 
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likelihood (£) of two hypotheses: (1) that the data were described by the standard LAT diffuse 
background models without any point sources (Ho), and (2) that the data were described by the 
existing model with an additional free parameter corresponding to the flux of a source at the target 
location (Hi). Utilizing the likelihood ratio test, we defined a test statistic (TS): 

Ts =- 2b (ffD- < n > 

After generating a map of the test statistic over the entire high-latitude sky, we iteratively 
refit regions around potential source candidates more carefully. The flux normalizations, spectral 
indices, and emission centroids of candidate sources with TS > 16 were refined while incorporating 
the flux normalizations of diffuse backgrounds and other candidate point sources with TS > 16 
within the ROI as free parameters in the fit. After refitting, only candidate sources with TS > 24 
were accepted into the list of source candidates 4 . Finally, to avoid duplicating sources in the 1FGL, 
we removed candidate sources with 68% localization errors overlapping the 95% error ellipse given 
for 1FGL sources. 

Our search of the high-latitude sky revealed 710 candidate sources, of which 154 were not in 
the 1FGL (36 of these candidate sources subsequently appear in the Second LAT Source Cata- 
log (Abdo et al. 2011)). We did not expect to recover all 806 high-latitude 1FGL sources, since 
the 1FGL is a union of four different detection methods and external seeds from the BZCAT and 
WMAP catalogs (Abdo et al. 2010b). However, since Sourcelike fits each energy bin indepen- 
dently, we expected to find source candidates that were excluded from the 1FGL, either because 
they had non-power-law spectra or they had hard spectra with too few photons to pass the 1FGL 
spectral analysis. We sacrificed some sample purity for detection efficiency in our candidate source 
list because stringent cuts on spatial extent and spectral shape were later applied. We obtained 
a final list of 385 high-latitude unassociated LAT sources and source candidates by combining the 
231 unassociated sources in the 1FGL with these 154 non-lFGL candidate sources. 

To check for consistency with the source analysis of the 1FGL, we performed an unbinned 
likelihood analysis with gtlike assuming that the unassociated sources were point-like with power- 
law' spectra. Our fitted fluxes and spectral indices are in good agreement wdth those in the 1FGL for 
the 231 unassociated 1FGL sources. The values are plotted in Figure 2, where it can be seen that 
the unassociated LAT sources span a wide range of fluxes and spectral indices. This wide range 
was taken into account when designing selection criteria for candidate dark matter satellites. The 
strong correlation between spectral index and flux is due to the improvement of the point-spread 
function (PSF) of the LAT with increasing energy and the relatively soft spectral dependence of the 
Galactic diffuse background. It is apparent that there are more non-lFGL source candidates in this 
sample with very hard spectra (spectral index ~ 1.0) and very low r fluxes (~ 10 -lu phcm~ 2 s _1 ). 
In the Appendix, we show’ that these source candidates are very likely spurious. 


4 Monte Carlo simulations have shown that 1 in 10 4 background fluctuations will be detected at. TS > 24 when fit 


with Sourcelike. (Wang 2011). 
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3.3. Spatial Extension Test 

The LAT has the potential to resolve some dark matter satellites as spatially extended 7 -ray 
sources. While the bulk of satellites are not spatially resolvable by the LAT, spatial extension is an 
important feature for distinguishing large or nearby satellites from point-like astrophysical sources 
(see Section 2). Assuming that the spatial and spectral distribution of 7 rays produced from dark 
matter annihilation factorize, the shape of the projected dark matter distribution (Equation (2)) 
is convolved with the LAT PSF. For an NFW dark matter distribution (Equation (5)) with scale 
radius r s at a distance D, the angular extent of a satellite can be characterized by the parameter 
a o — r s /D. Approximately 90% of the integrated J-factor comes from within the angular radius 
ao (Strigari et al. 2007) and the LAT is sensitive to the spatial extension of satellites with a > 0.5 °. 

We used the likelihood ratio test, as implemented by Sourcelike, to test sources for spatial 
extension. We defined a test statistic for extension as 

TS„ = -2i„(lAe=d) (12) 

VM-nNFWj / 

= TSnFW - TSpoint (13) 

where TS p0 mt was the test statistic of the candidate source assuming that it had negligible extension 
(ao much smaller than the LAT PSF) and TSnfw was the test statistic of the candidate source 
when ao was fit as a free parameter. In both cases, the position of the source was optimized during 
the fit. 

We sought to define a cut on the value of TS ext to eliminate 99% of point sources over the 
range of spectral indices and fluxes found in the unassociated LAT sources. This cut was labeled 
TSg xt - While the point and extended hypotheses are nested and TS ext is cast as a likelihood ratio 
test, it is unclear whether this analysis satisfied all of the suitable conditions for the application the 
theorems of Wilks (W r ilks 1938) or Chernoff (Chernoff 1954). Therefore, we relied on simulations 
to parameterize TSg Xt as a function of source flux and spectral index. 

To evaluate TSf xt over the pertinent range of source fluxes and spectral indices, the unassoci- 
ated LAT sources were bracketed with 10 representative power-law models (Table 1 and blue stars 
in Figure 2). For each of the 10 representative models. 1000 independent sources were simulated 
at random locations in the high-latitude sky using the LAT simulation tool, gtobsslm, and the 
spacecraft pointing history for our one- year data set. To accurately incorporate imperfect modeling 
of background point and diffuse sources, we embedded the simulated point sources in one year of 
LAT data and calculated TS ext for each simulated source. We defined TSg Xt for each representative 
model as the smallest value of TS ext that was larger than that calculated for 99% of simulated 
point sources. The values of TSg Xt for all 10 models were calculated independently of the TS > 24 
detection cut (Table 1). 

We used a bilinear interpolation to estimate the value of TSg Xf for any point in the space 
spanned by the grid of flux and spectral index. Since each measurement of source flux and spectral 



13 - 


index has a statistical uncertainty, we interpolated to the largest value of TS^t that was consistent 
with the ±lcr error for each source for a conservative estimate of TSg® t . 

3.4. Spectral Test 

To select sources that were spectrally consistent with WIMP dark matter annihilation, we 
designed a test for spectral curvature. The continuum 7 -ray emission from WIMP annihilation 
has two different contributions: secondary photons from tree-level annihilation (Baltz et al. 2008; 
Cesarini et al. 2004) and additional photons from QED corrections - i.e., final state radiation 
(FSR) (Beacom et al. 2005). For tree-level annihilations, the leading channels among the kine- 
matically allowed final states are predicted to be bb, it, W + W~~ , Z°Z°, and t + t~. The 7 -ray 
spectra from these channels are quite similar, except for the r-channel which is considerably 
harder (Cesarini et al. 2004). We chose the bb channel as a representative proxy for the tree-level 
annihilation spectrum. 

We defined a test statistic to evaluate the consistency of the data with dark matter and power- 
law spectra. This spectral test statistic, 

TS spec — ~2\n 

= TS^ — TSp W i (14) 

was the difference in source TS calculated with an unbinned analysis using gtlike assuming a bb 
dark matter spectral model (TS^j) and a power-law (TS pw i) spectral model. These two hypotheses 
are not nested, and thus the significance of this test was evaluated with simulations. When per- 
forming our fits, we modeled candidate sources as point-like 5 and left their fluxes and the fluxes of 
the diffuse backgrounds free. Additionally, the power-law and dark matter spectral models contain 
a spectral free parameter (the dark matter mass or power-law index). 

Using the same representative simulations described in Section 3.3, we defined TS® pec to be 
the value of TS spec which was larger than that calculated for 99% of simulated power-law sources 
(Table 1). When calculating TSg pec for a particular source, we chose the largest value from a 
bilinear interpolation to the ±ler errors on fitted flux and spectral index (as discussed at the end of 
Section 3.3). These tests of spatial extension and spectral character allowed us to select non-point- 
like and non-power-law sources with a contamination of 1 in 10 4 assuming they were independent. 

‘'From simulations, this has been found to be a conservative way to estimate TS spec (Wang 2011). 
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4. Results 

4.1. Search for Dark Matter Satellites 

We applied the cuts on spatial extension and spectral character to select DM satellite candi- 
dates from the 385 unassociated high-latitude LAT sources and source candidates. Two of the 385 
unassociated sources, 1FGL J1302.3— 3255 and 1FGL J2325.8— 4043, passed the cut on spatial ex- 
tension. One of these, 1FGL J1302.3— 3255, also passed our spectral test, preferring a bb spectrum 
to a power-law spectrum. However, we do not believe that either of these sources is a viable DM 
satellite candidate for reasons discussed below. 

As their names imply, both 1FGL J1302.3— 3255 and 1FGL J2325.8— 4043 are present in the 
1FGL (Abdo et al. 2010b); we summarize information about each source in Table 2. While 
1FGL J1302. 3—3255 was unassociated when the 1FGL was published, and has previously been con- 
sidered as a promising dark matter satellite candidate (Buckley & Hooper 2010), it has since been 
associated with a millisecond pulsar by radio follow-up observation (Hessels et al. 2011). The other 
possibly extended source, 1FGL J2325.8— 4043, has a high probability of association with two AGN 
in the first LAT AGN Catalog (Abdo et al. 2010f). 1FGL J2325.8— 4043 is assigned a 70% prob- 
ability of association to 1ES 2322— 409 and a 55% probability of association with PKS 2322— 411. 
Additionally, each of these sources was checked for consistency with an E~ x power-law spectrum, as 
would be expected from some FSR models (Essig et al. 2009), and both can be excluded with 99% 
confidence (Wang 2011). We further discuss sources with hard spectral indices in the Appendix. 

Since AGN are not expected to be spatially extended at an angular scale resolvable by the LAT, 
we cross checked 1FGL J2325.8— 4043 against the Second LAT Source Catalog (Abdo et al. 2011). 
In two years of data, two sources were found within 0.5 0 of the location of 1FGL J2325.8— 4043. In 
one year of data, these two sources could not be spatially resolved, but their existence was enough 
to favor an extended source hypothesis. Spurious measurements of a finite extent for point sources 
are not unexpected. Testing extension with a purity of 99%, the Poisson probability of finding at 
least one spurious source in our 385 tests is 98%. Since 1FGL J1302.3— 3255 is associated with a 
pulsar and 1FGL J2325.8— 4043 does not appear to be truly extended, we conclude that there w'ere 
no unassociated, high-latitude spatially extended 7-ray sources in the first year of LAT data. Thus, 
according to the criteria defined in Section 3, no viable dark matter satellite candidates were found. 


4.2. Contamination from Pulsars 

To better understand possible misidentification of pulsars (such as 1FGL J1302.3— 3255) as 
candidate dark matter satellites, we applied our test of spectral shape to 25 high-latitude LAT 
detected pulsars. Of these 25 pulsars, 14 were identified when the 1FGL was published and 11 were 
subsequently identified through follow-up observations by radio telescope (Hessels et al. 2011). 

Interestingly, 24 of these pulsars passed our spectral cut, preferring a bb spectrum to a power- 
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law spectrum. This can be understood by comparing the exponentially cutoff power-law model 
commonly used to fit pulsars with a bb WIMP annihilation spectrum. The exponentially cutoff 
power law has the form (Abdo et al. 2010g): 


dN rr-r , E 
-=K£ &v ex p(- — 


(15) 


where T is the photon index at low energy, E cut is the cutoff energy, and If is a normalization 
factor (in units of phcm - 2 s -1 MeV -1 ). In Figure 3, we plot both the exponentially cutoff power- 
law model and a low-mass (Mwimp ~ 25GeV) bb spectrum and show that for E > 200 MeV the 
two curves are very similar. 

By fitting bb spectra to the 25 LAT-detected pulsars, we found that they tend to be best 
fit by low dark matter masses (Figure 4). Although our statistics are limited, the distribution 
peaks around 30GeV, with most pulsars having a best- fit dark matter mass Mwimp < 60GeV. 
This suggests that unidentified, high-latitude pulsars can present a source of confusion in spectral 
searches for dark matter satellites. In general, many unassociated LAT sources have spectra that 
are inconsistent with a power-law model (Bonamente 2010; Abdo et al. 2011). The fact that these 
sources passed our spectral test does not imply that they are best fit by bb spectra, merely that bb 
spectra fit better than a simple power law. These unassociated, non-power-law sources were not 
found to share a consistent spectrum, as would be expected from dark matter annihilation. 


The abundance of non-power-law 7 -ray sources emphasizes the importance of testing for spatial 
extension when attempting to identify dark matter satellites at high latitudes. Some concerns 
remain due to the fact that the LAT detects spatially extended pulsar wind nebulae located around 
some pulsars (Ackermann et al. 2011b). However, we do not expect the older pulsars at high 
Galactic latitudes to have nebulae that are spatially extended on a scale detectable by the LAT. 
Of course, there is always a risk that a chance coincidence with a low-flux neighboring source will 
cause apparent source extension. 


5. Interpretation in the Context of A r -body Simulations 

No high- latitude unassociated LAT source candidates passed our dark matter satellite selection 
criteria. This is combined with the simulations in Section 2 to constrain a conventional 100 GeV 
WIMP annihilating through the bb channel. Monte Carlo simulations were used to determine the 
detection efficiency for dark matter satellites as a function of flux and spatial extension. For multiple 
realizations of each IV-body simulation, we calculated the probability of detecting no satellites from 
the detection efficiency of each simulated satellite. Averaging over these simulations and increasing 
( 1 av ) until the probability of detecting no satellites drops below 5%, we were able to set a 95% 
confidence upper limit on (crv). 
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5.1. Detection Efficiency 

The detection efficiency of our selection was defined as the fraction of dark matter satellites 
that pass the cuts in Section 3 and was calculated from Monte Carlo simulations. The efficiency for 
detecting a dark matter satellite depended on spectral shape (i.e., dark matter mass and annihilation 
channel), flux, and spatial extension. For a 100 GeV WIMP annihilating through the bb channel, we 
examined the efficiency for satellites with characteristic fluxes ranging from 5.0 x 10” 11 phcm~~ 2 s" 1 
to 5.0 x 10~ 8 ph cm~ 2 s^ 1 and characteristic spatial extension (as described in Section 3.3) from 
0.5° to 2.0°. These ranges were chosen to reflect the fluxes of the unassociated high-latitude LAT 
sources and angular extents to which the LAT is sensitive 6 . 

For each set of characteristics listed in Table 3, we utilized gtobssim to simulate 200 dark 
matter with NFW profiles and bb spectra from a 100 GeV WIMP. These simulations were embedded 
in LAT data at random high-latitude locations, and Sourcelike was used to compute TS ext , TS spec , 
and the detection TS for each. The satellite detection efficiency was computed as the fraction of 
satellites with Sourcelike TS > 24, TS ext > TSg Xt , and TS spec > TSg pec . The first requirement 
was included as a proxy for the efficiency of the source finding algorithm. The creation of this 
efficiency table (Table 3) was computationally intensive and the result is likely model dependent, 
which limited this analysis to the examination of only the 100 GeV bb model. To expedite the 
generation of this table, we found the flux value with efficiency < 0.05 and conservatively set the 
efficiency for sources with less flux to 0. 


5.2. Simulated Satellite Distributions 

The VL-II and Aquarius simulations (described in Section 2) were used to predict the Galactic 
dark matter satellite population. Picking a vantage point 8.5 kpc from the center of each simulation 
(the solar radius), we calculated the spatial extension and integrated /-factor for each satellite. To 
account for variation in the local satellite population, we repeated this procedure, creating six 
realizations from maximally separated vantage points at the solar radius. It is important to note 
that while the VL-II and six Aquarius simulations are statistically independent, these different 
realizations are not - i.e., the same satellites appear in multiple realizations. Thus, we have seven 
independent simulations, each with six not-independent realizations (collectively referred to as 42 
“visualizations” ) . 

After excluding undetectably faint satellites with /-factors an order of magnitude less than the 
lower bound the /-factor of Draco (Section 2.3), we compared the distributions of /-factors and 
spatial extensions across the visualizations. We found reasonable agreement in the detectable satel- 
lite distributions between the VL-II and Aquarius simulations. The variation in satellite number 


“The 68% containment radius of the LAT PSF, which depends on photon energy and angle of incidence, can be 
approximated by the function, 0.8°(nJ|v) ' * (Abdo et al. 2010b), yielding ~ 0.8° at lGeV and ~ 0.13° at 10 GeV. 
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in bins of flux and spatial extension was much larger between different simulations than between 
realizations of the same simulation (as is expected, because the realizations are not independent). 
On average, satellites with oo > 0.5 ° make up ~ 30% of the total integrated J-factor from satellites 
in these simulations. 


5.3. Upper Limits 

For each of the 42 visualizations of VL-II and Aquarius, we calculated the 7 -ray fluxes of all 
satellites for a given ( av ) using Equation (1). With these fluxes and the true spatial extension for 
each satellite, we performed a bilinear interpolation on Table 3 to determine the detection efficiency 
for each satellite. The probability that the LAT would observe none of the satellites in visualization 
i is 

p i{{w )) = n ^ 1 - tiAi™))) ( ie ) 

3 

where €ij is the detection efficiency for satellite j in visualization i. Because there is no reason to 
favor any one visualization, we calculated the average null detection probability over the N — 42 
visualizations as 

1 N 

( 6 * * * * * * * * * * 17 ) 

i 

To set an upper limit on the dark matter annihilation cross section, we increased (av) until the 
probability of a null observation was < 5%, i.e. P < 0.05. This corresponds to 95% probability that, 
for this (av ) , at least one satellite would have passed our selection criteria. Using this methodology, 
the LAT null detection constrains (av) to be less than 1.95 x 10~ 24 cm 3 s _1 for a 100 GeV WIMP 
annihilating through the bb channel. 


6. Discussion and Conclusions 

We performed a search for dark matter satellites using one year of Fermi -LAT data. After 

completing an independent search for 7 -ray sources, we constructed tests to evaluate both source 

extension and spectral shape. We distinguished dark matter satellite candidates by selecting spa- 

tially extended sources with 7 -ray spectra consistent with those produced by dark matter particle 

annihilation to bb. Our initial scans selected two potentially extended sources; however, follow-up 

analyses revealed that neither of them is a valid dark matter satellite candidate. Therefore, we 

concluded that, given our pre-defmed search criteria, there were no signals of dark matter satellites 

in the first year LAT data. 

Using ACDM-based theoretical predictions from the Aquarius and Via Lactea II numerical 

simulations of the Galactic dark matter distribution, we estimated the number of dark matter 

satellites that could be observed assuming a 100 GeV bb model with varying cross section. We 
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quantified the detection efficiency for these satellites, and used it to set an upper limit on the velocity 
averaged annihilation cross section of 1.95 x 10~ 24 cm 3 s“ 1 for a 100 GeV WIMP annihilating to 
bb. This limit is approximately 60 times greater than the expected value of the thermal relic cross 
section, and it is about an order of magnitude less stringent than the best one-year limit from known 
dwarf spheroidals (Abdo et al. 2010e). This difference in upper limits can be accounted for by the 
fact that the analysis of dwarf spheroidals effectively probes the flux limit of the LAT detector, 
while our selection on spatial extension and spectral shape limits us to higher source fluxes. 

We have presented a novel technique to search for dark matter satellites using the LAT. Further 
data will improve the prospects for finding dark matter satellites using the methodology developed 
in this paper. Higher quality data and analysis techniques are likely to reveal the presence of more 
unassociated sources. Improved analysis techniques will make it possible to select dark matter 
candidates with increased efficiency while maintaining the same discrimination power. Addition- 
ally, using the analysis in this paper as a guide, one can consider a broader class of dark matter 
annihilation models. 
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A. Investigating High Latitude Spurious Sources 

In this Appendix, we estimate the number and distribution of spurious sources arising from 
our source search method. The TS of each source can be related to the probability that such an 
excess can be obtained from background fluctuations alone. It was shown (Mattox et al. 1996) 
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that, for the EGRET 7 -ray telescope, the distribution of TS at a fixed point for background-only 
simulated data followed a y 2 / 2 distribution with one degree of freedom as expected from Chernoff’s 
theorem (Chernoff 1954). Once convinced that the necessary conditions were met for Chernoff’s 
theorem, it was possible to use standard combinatorial arguments, given the method used to scan 
the sky, to calculate the false source rate at any TS value. However, in our analysis it was unclear 
whether all of the suitable conditions of Chernoff’s theorem were satisfied. Indeed, Wang (2011) 
shows that, in our case, the source TS distribution does not follow that predicted by Chernoff’s 
theorem. 

To estimate the number of spurious sources in our candidate source list, we applied our source 
search procedure to a one-year background-only Monte Carlo simulation (Galactic and isotropic 
diffuse only). The one-year LAT observation simulation was generated using gtobssim with the 
actual LAT pointing history and application of the same event selection (P6_V3„DIFFUSE IRFs, 
time range, energy range, zenith angle and rocking angle). 

We searched the background-only simulated data at \b\ >20° using Sourcelike, and found 
193 candidate sources with TS > 24. Of course, in this case all 193 sources were spurious. This 
distribution was uniform over the high-latitude sky and was consistent with arising from random 
fluctuations in the background. The locations of the spurious sources have no coincidences with 
any of the 385 unassociated sources and source candidates found in the LAT data (see Section 3.2). 

The 193 spurious sources were individually analyzed with gtlike, using a point-source spatial 
model and a power- law spectral model. Figure 5 shows the distribution of the spectral index and 
the integral flux from 200 MeV to 300 GeV for these spurious sources (blue crosses) and the 385 
unassociated sources and source candidates (black and red squares). The same correlation between 
spectral index and flux holds for both the Monte Carlo simulation and LAT data. 

We specifically studied the spurious sources and unassociated sources with spectral indices in 
the range from 1.0 to 0.7 since a cluster of sources in this range are predicted by some leptophilic 
dark matter models (Arvanitaki et al. 2009; Grasso et al. 2009). In this range of spectral indices, 
there are 12 spurious sources in the Monte Carlo simulation and 13 unassociated source candidates 
in the LAT data. We applied a more stringent set of photon selection cuts, to obtain a cleaner 
photon sample than the “Diffuse” class (Abdo et al. 2010c). This cut removed most of the residual 
high-energy charged-particle background from the regions of the 13 unassociated source candidates 
and resulted in a general reduction of source significance and the disappearance of two sources 
(TS ~ 0). This implied that the detection of some of these unassociated hard source candidates 
was due to high-energy charged particle contamination. There were no obvious differences in the 
spatial or spectral distribution of stringently selected photons in the regions of the 13 unassociated 
source candidates and simulated photons from the 12 spurious sources. Additionally, we reanalyze 
these 13 unassociated source candidates using 25 months of LAT data, and found that 10 of the 
13 sources decrease in significance and none of them increased substantially. Given the discussion 
above, the 13 unassociated source candidates with spectral indices in the range from 1.0 to 0.7 are 



consistent with being spurious and do not provide evidence for dark matter satellites. 
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Mass (M 0 ) 

Fig. 1. — Distribution of satellite mass and distance for the original VL-II satellites (in black) and 
the extrapolation to low-mass satellites (in red). Lower J-factors reside in the upper left while 
higher J-factors lie to the lower right. Contours of constant J-factor (J ex ^-) ran from the 
upper right to the lower left. One such contour is shown for the Draco dwarf spheroidal galaxy 
assuming a mass of 10 8 Mq at a distance of SOkpc. Satellites lying in the hatched region above this 
line have lower J-factors than that of Draco. 
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Fig. 2. — Distribution of spectral indices and integral fluxes from 200 Me V to 300 GeV for the 385 
high-latitude unassociated sources and source candidates. The squares are the 231 unassociated 
sources from the 1FGL catalog, while the triangles are the 154 additional source candidates detected 
with Sourcelike. The stars are the 10 representative porver-law models in Table 1. 
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Fig. 5. — Distribution of spectral indices and integral fluxes from 200 MeV to 300 GeV for the 385 
high-latitude unassociated sources and source candidates. The squares are the 231 unassociated 
sources from the 1FGL catalog, while the triangles are the 154 additional source candidates detected 
with Sourcelike. The circles are the 193 spurious sources found in a Monte Carlo simulation of 
background only. 



Table 1. Values for TS*® t and TSf® ec 


Model Number 

Spectral Index 

Flux'"' 

( phcm" 2 s -1 ) 

CPq 99 
1 ^ext 

rpq99 
L °spec 

1 

0.9 

2.0 x 10~ 10 

6.18 

2.38 

2 

0.9 

8.0 x 10~ n 

7.87 

2.46 

3 

1.5 

1.1 x 10~ 9 

5.09 

4.96 

4 

1.5 

2.0 x 10~ 10 

14.98 

2.88 

5 

2.0 

1.2 x 10~~ 8 

5.11 

2.24 

6 

2.0 

1.2 x 10~~ 9 

9.63 

4.28 

7 

2.5 

2.1 x 10 ~ 8 

6.74 

1.78 

8 

2.5 

0.5 x 10~ 8 

10.78 

5.66 

9 

3.0 

1.7 x 10~ 8 

9.81 

2.14 

10 

3.0 

1.0 x 10~ 8 

11.87 

6.02 


'“'integral flux from 200 MeV to 300 GeV. 


Note. — Cut values excluding point-like sources (TS 99 t ) and sources 
with power-law spectra (TS 99 e c) at 99% confidence. These values are 
independent of the TS > 24 cut and were calculated from Monte Carlo 
simulations of the 10 typical power-law point source models. 


Table 2. Two Candidate Extended Sources 


Source ID 

l, b M 
(deg) 

Flux' 6 ' 

(phcm _2 s _1 ) 

ao (c) 

(deg) 

TS ext 

TSspcc 

1FGL J1302.3— 3255 

305.58, 29.90 

1.33 x 10~ 8 

1.2 

9.3 

4.6 

1FGL J2325.8— 4043 

349.83, -67.74 

2.12 x 10 -8 

1.3 

13.2 

-19.1 


'“'Best-fit source position (Abdo et al. 2010b). 

^'integral flux from 200 MeV to 300 GeV interpolated from best-fit power-law 
model (Abdo et al. 2010b). 

'“'Best-fit spatial extension and ±la error using an NFW profile. 
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Table 3. Detection Efficiency 


Flux (a) 

(ph cm” 2 s -1 ) 

0.5° 

Extension 

1.0° 

2.0° 

0.2 x 10” 8 

< 0.05 

< 0.05 

< 0.05 

0.5 x 10” 8 

0.16 

0.28 

0.31 

1.0 x 10” s 

0.74 

0.76 

0.83 

2.0 x 10” 8 

0,99 

1.0 

0.99 

5.0 x 10“ s 

1.0 

1.0 

1.0 


^Integral flux from 200 MeV to 300 GeV. 

Note. — Satellite detection efficiency for 100 GeV 
WIMP annihilating through the 66 channel. 



